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We report the results of a search for pair production of doubly-charged Higgs bosons via pp — > 
H ++ H~~ X — > fi + pi + fi~ fi~ X at y/s — 1.96 TeV. We use a dataset corresponding to an integrated 
luminosity of 1.1 fb _1 collected from 2002 to 2006 by the DO detector at the Fermilab Tevatron 
Collider. In the absence of an excess above the standard model background, lower mass limits 
of M(J/± ± ) > 150 GeV/c 2 and M(ff± ± ) > 127 GeV/c 2 at 95% C.L. are set, respectively, for left- 
handed and right-handed doubly-charged Higgs bosons assuming a 100% branching ratio into muons. 
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PACS numbers: 14.80.Cp, 13.85.Rm 



In the standard model (SM) of electroweak interac- 
tions, elementary fermions and bosons acquire mass via 
a weak isospin scalar doublet. This mechanism results 
in the existence of an additional particle, the Higgs bo- 
son, which has not yet been observed. Extensions of the 
Higgs sector involving higher isospin multiplets predict 
the existence of doubly-charged Higgs bosons which can 
be relatively light and hence accessible at current ex- 
perimental facilities. Doubly-charged Higgs bosons ap- 
pear in many scenarios such as left-right symmetric mod- 
els fil, Higgs triplet models @, and Little Higgs mod- 
els [3j. At the Fermilab Tevatron Collider, the two 
main production mechanisms are pair production via 
pp — > Zj^*X — > H ++ H X and single production via 
WW fusion, pp -► W ± W ± X -► H ±:t X. However, 
higher isospin Higgs multiplets are generally severely con- 
strained by p = m^r / (cos 6 w m z) 2 = 1 at tree level. The 
existing phenomenological and theoretical constraints are 
easily satisfied when the W^W^ — » i? ±± coupling is 
vanishing [4j. If the H ++ coupling to W boson pairs 
is suppressed, the dominant final states are expected to 
be like-sign lepton pairs. Left-handed (Hf~^) and right- 
handed (Hft ) states are distinguished by their coupling 
to left-handed and right-handed leptons, respectively. 
The pair production cross section for left-handed doubly- 
charged Higgs bosons for 100 < M(i? ±:t ) < 200 GeV/c 2 
is about a factor two larger than that for the right-handed 
states due to different couplings to the intermediate Z bo- 
son Previous searches for iJ ±± have been performed 
by the LEP collaborations 0] in e + e~ collisions and by 
the DO [7] and CDF [|[ collaborations at the Tevatron 
pp collider. This Letter presents the results of a direct 
search for pp — ► H ++ H X with i? ±± — > 
DO collaboration with improved sensitivity. 



H fi^ by the 



The main DO detector systems are a central tracking 
system, a liquid-argon and uranium calorimeter, and a 
muon detector [9(. The central tracking system consists 
of the silicon microstrip tracker (SMT) and the central 
fiber tracker (CFT) surrounded by a 2 T solenoidal mag- 
net, with designs optimized for tracking and vertexing 
capability at pseudorapidity [nj \tj\ < 3 and \rj\ < 2.5, 
respectively. The liquid-argon and uranium calorimeter 
has a central calorimeter (CC) covering a region up to 
|ry| « 1.1 and two end calorimeters (EC) extending the 
coverage to \r]\ sa 4.2, with each housed in a separate 
cryostat llj. The muon detector has layers of propor- 
tional drift tubes and scintillation counters before and 
after a 1.8 T iron toroid [12]. This analysis is based on 
the Run II data set collected with the DO detector at 
the Fermilab Tevatron Collider at y/s — 1.96 TeV from 
April 2002 to February 2006 corresponding to 1.1 fb _1 . 
Events are collected using a suite of dimuon and single 
muon triggers. 



In the previous DO analysis Q, two like-sign muons 
were required in the final state. In this analysis, we re- 
quire a third muon, which increases the sensitivity by 
decreasing backgrounds. We follow five steps to select 
events. In the first step (SI), events are required to have 
at least two muons. Each muon must have a transverse 
momentum px > 15 GeV/c and \r]\ < 2.0. Muons are 
selected using patterns of hits in the wire chambers and 
scintillators in the muon system. Each muon must be 
matched to a track in the central tracker with at least 
five hits in the CFT layers and at least two hits in the 
SMT layers. Muons from cosmic rays are removed by 
using a timing information on the hits in the scintillator 
layers. 

In the second step (S2), isolation criteria based on 
the calorimeter and tracking information are applied 
to remove the background from multijet production 
with muons originating from in-flight decay of pions or 
kaons, or from semi-leptonic decays of B or D mesons. 
The sum of the transverse energies of the calorime- 
ter cells in an annulus of radius 0.1 < 7Z < 0.4, where 
1Z = W (A0) 2 + (Ary) 2 and (j> is the azimuthal angle, 
around the muon direction is required to be less than 
2.5 GeV. A similar condition is defined for the scalar 
sum of the pr of all tracks, excluding the muon in a cone 
of radius 1Z = 0.5 centered around the muon, which must 
be less than 2.5 GeV/c. 

Selection S3 reduces the remaining Z — > and 
multijet backgrounds. The azimuthal angle A(f> between 
at least one pair of muons is required to be less than 
2.5 radians, since the two muons from Z boson decays 
are mostly back-to-back. This requirement also rejects a 
fraction of the multijet background with nearly back-to- 
back muons. 

Selection S4 requires at least two muons to be of like 
sign. The final selection (S5) requires a third muon, sat- 
isfying the SI selection and the isolation selection criteria 
S2 but without the minimum hit requirement on the cen- 
tral track. 

The dominant background in this analysis arises from 
electroweak processes where real high pt muons are cre- 
ated from W or Z boson decays as well as non-isolated 
muons originating from jets. The SM backgrounds and 
signal processes are generated with pythia [l]| and 
normalized using the theoretical cross section. The 
Z/7* — * l + t~ cross section is calculated at next-to- next- 



to-leading order (NNLO) [14j. The tt cross section is cal- 
culated at NNLO [H and the WW.ZZ and WZ cross 
sections are calculated with MCFM [l6| at next-to-leading 
order (NLO). All samples are processed through the DO 
detector simulation based on GEANT 17] and the same 
reconstruction software as for the data. The muon recon- 
struction and isolation efficiencies differ between Monte 
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TABLE I: The expected numbers of events for a signal with M(i2J~ ) = 140 GeV/c 2 and background and the number of 
observed events after each selection step. The statistical and systematical uncertainties are combined in the table. 



Selection 


Preselection 


Isolation 


Ac/> < 2.5 


Like sign 


Third muon 




SI 


S2 


S3 


S4 


S5 




69181 ±4642 


58264 ±3910 


4936 ± 333 


5.3 ±1.6 


< 0.01 


Multijet 


4492 ±120 


194 ±18 


18-1-9 
lo zrz z 


D.o zrz U.o 


9-1-01 
U.z zrz u.l 




328 ± 25 


269 ±21 


20 ±3 


< 0.01 


< 0.01 


ti 


38 ±3 


20 ±1 


14 ±1 


0.03 ±0.01 


< 0.01 


WW 


40 ±3 


34 ±2 


20 ±1 


< 0.01 


< 0.01 


wz 


19 ±1 


16 ±1 


11 ± 1 


2.95 ±0.20 


1.62±0.11 


zz 


10 ±1 


9±1 


5±1 


0.63 ±0.05 


0.47 ±0.03 


Total background 74108 ± 4644 58806 ± 3910 


5024 ± 333 


15.2 ±1.8 


2.3 ±0.2 


Signal 


20.5 ±2.7 


18.5 ±2.4 


16.3 ±2.1 


11.6 ±1.5 


10.1 ±1.3 


Data 


72974 


58763 


4558 


16 
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Carlo (MC) and data, and these differences are corrected. 
Trigger efficiency corrections are not applied to the MC 
sample. Instead, the MC sample are normalized to the 
data using the Z boson mass peak at the selection level 
S2. 

Another important background comes from multijet 
production, mainly bb events decaying semi-lcptonically 
into muons that appear isolated. The multijet back- 
ground is derived from the data sample with non-isolated 
muons obtained by inverting the isolation requirements 
for both muons after the selection SI. The efficiency 
of the isolation requirement is assumed to be identical 
for multijet events with like-sign and opposite-sign muon 
pairs. It is also assumed that all like-sign events after 
subtracting SM backgrounds are multijet events. The 
SM backgrounds are subtracted in the following samples 
used for the multijet background determination. The 
total number of multijet events before the isolation re- 
quirement (4492 ± 120) is then given by the number of 
non-isolated events for all charge combinations multiplied 
by the ratio of the total number of events to the num- 
ber of non-isolated events in the like-sign sample. The 
number of multijet events after the isolation requirement 
(194 ± 18) is obtained by multiplying this number with 
the isolation efficiency (4.3 ±0.5)%, given by the ratio of 
isolated to all like-sign multijet events. 

A second instrumental background arises from Z/"f* — > 
events in which the charge of one of the muons is 
misidentified. The first source of charge misidentifica- 
tion is due to fewer CFT layers at large r\ and a conse- 
quent increase in the charge misidentification probabil- 
ity. The second source affects very high pr tracks for 
which the uncertainty on the measured curvature can 
cause charge misidentification. The charge misidentifica- 
tion rate is obtained by dividing the number of like-sign 
events (SI, S2 and S4) by the number of events with- 
out the like-sign requirement (SI and S2) in the dimuon 
invariant mass region above 70 GeV/c 2 , after subtract- 
ing the SM sources of background except Z/j* — > fi + n~ 



events from the data. This mass requirement removes 
most multijet background events in the low mass range. 
From these ratios, we determine the average probabil- 
ity for charge misidentification in data and MC to be 
P data = (6.2 ± 1.1) x 10- 4 and P MC = (3.1 ±0.4) x 10" 4 , 
respectively, assuming the multijet background is negli- 
gible. The uncertainties are statistical. Since the charge 
misidentification rate in MC is underestimated, the ratio 
of -Pdata to Pmc is taken as a correction equal to 2.0±0.4. 
This ratio is applied to the Z/j* — > [i + f-i~ MC sample 
when estimating the like-sign contribution. 

The distributions of dimuon invariant mass and A</> 
after the selection SI are shown in Fig. [1] (a) and (b). 
The data are compared with the sum of the background 
contributions. For those events with more than one pair 
of muons fulfilling the selection criteria, the dimuon in- 
variant mass and A<j) are calculated only for the pair 
with the highest individual momenta. The numbers of 
remaining events after each selection are shown in Ta- 
ble |U There is good agreement between data and the 
sum of the backgrounds. Figure Q] (c) and (d) show the 



TABLE II: The numbers of observed and expected back- 
ground events after each selection criterion with the like-sign 
requirement applied together with SI. The statistical and 
systematical uncertainties are combined in the table. 



Selection 


Preselection 


Isolation 


Ac/> < 2.5 


(Like-sign) 


SI & S4 


S2 


S3 


Z/ 7 * -> 


84 ±24 


42 ±12 


5.3±1.6 


Multijet 


1620 ± 34 


70 ±5 


6.3±0.8 


Z/j* -» T+T~ 


3.2 ±1.3 


0.2 ±0.3 


< 0.01 


tt 


6.6 ±0.5 


< 0.1 


< 0.1 


WW 


0.08 ±0.02 


0.04 ±0.01 


< 0.01 


wz 


5.14 ±0.35 


4.25 ±0.29 


2.95±0.20 


zz 


1.12 ±0.08 


0.90 ±0.06 


0.63±0.05 


Total background 


1720 ±41 


117±13 


15.2±1.8 


Data 


1678 


96 


16 
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FIG. 1: Distributions of the dimuon invariant mass and A<f> 
between the two muons for data compared to the sum of MC 
backgrounds after the selection SI (a, b), the preselection 
SI with the like-sign requirement S4 (c, d) and the final re- 
quirement S1-S5 (e, f). The signal expected for a left-handed 
H ±zk , with M(H ±:k ) = 140 GeV/c 2 , is also shown by the 
open histogram (e, f). 



dimuon invariant mass and A<f> distributions after the SI 
and S4 requirements. The excess of events at 150 GeV/c 2 
has a significance of less than 2.6a. Table HT1 gives the in- 
dividual like-sign backgrounds after the various selection 
stages. This demonstrates that the like-sign backgrounds 
are well understood. 

After all five selection criteria, three data events re- 
main, in good agreement with the SM background ex- 
pectation of 2.3 ± 0.2 events. Total signal efficiencies 
are 32%-34% and are nearly independent of mass. The 
dimuon invariant mass and A</> distributions for these 
events are compared to the sum of the backgrounds in 
Fig. [T](e) and (f). 

Since no excess is observed, we use the dimuon invari- 
ant mass distribution in Fig.[T](e) to compute upper lim- 
its on the production cross section times branching frac- 
tion as a function of M (iJ ±± ) using the CLs method (lij 
as implemented in the MCLIMIT program [19(. The ex- 
pected rate for the signal as a function of M(H ±:t ) is 
determined by the next-to-leading order (NLO) cross sec- 




D0, L = 1.1 fb" 1 

— Theory(L) 
-Theory(R) 

— Observed 
— - Expected 

±1 a band 



80 100 120 140 160 180 200 

(GeV/c 2 ) 

FIG. 2: The cross section limit as a function of the Higgs 
mass M(H ±:h ) at the 95% C.L. The mass regions excluded 
by LEP and CDF are also shown. The ±lcr uncertainty on 
the expected limit is given by the yellow band (color online). 



tion [B[ and measured luminosity, corrected for the signal 
efficiency. 

A number of systematic uncertainties on signal and 
background are taken into account in the limit calcula- 
tion. The uncertainties on the correction of the muon 
identification are 2% and 6% for backgrounds and signal, 
respectively. The uncertainty on the isolation efficiency 
for the multijet background is 12%. The 20% uncertainty 
on the correction for charge misidentification is included. 
The uncertainty on the luminosity for signal is estimated 
to be 6.1% [2(|. The uncertainty on the normalization 
using NNLO MC SM background production cross sec- 
tions is taken to be 5%. The PDF uncertainties on the 
cross section for backgrounds are taken to be 4% [2lj ]. 

The cross section limit as a function of M(H ±:t ) is 
shown in Fig.[2]together with the theoretical cross section 
for left- and right-handed doubly charged Higgs bosons. 
At the 95% C.L., lower mass limits of 150 GeV/c 2 for left- 
handed and 127 GeV/c 2 for right-handed doubly-charged 
Higgs bosons are obtained. This significantly extends the 
previous mass limit [§] for a doubly-charged Higgs boson 
decaying into muons. 
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